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Abstract 



Diluted arrays of many optical apertures will be able to provide high- 
resolution snapshot images if the beams are combined according to the 
densified-pupil scheme. We show that the same principle can also provide 
coronagraphic images, for detecting faint sources near a bright unresolved 
one. Recent refinements of coronagraphic techniques, i.e. the use of a 
phase mask, active apodization and dark-speckle analysis, are also ap- 
plicable for enhanced contrast. Implemented in the form of a proposed 
50-500m Exo-Earth Discoverer array in space, the principle can serve to 
detect Earth-like exo-planets in the infra-red. It can also provide images 
of faint nebulosity near stars, active galactic nuclei and quasars. Calcu- 
lations indicate that exo-planets are detectable amidst the zodiacal and 
exo-zodiacal emission faster than with a Bracewell array of equivalent 
area, a consequence of the spatial selectivity in the image. 
Keywords: extrasolar planets; instrumentation; image processing. 
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1 Introduction 



An early proposal for searching exo-planets involved coronagraphy aboard the 
"Large Space Telescope" (Bonneau et al. 1975) which became the Hubble Space 
Telescope (HST). It was shown that telescope rotation and synchronous image 
subtraction could remove much of the coronagraphic residue, a speckled halo of 
star light expected to outshine the planet image by 3 or 4 orders of magnitude. 
The corresponding Faint Object Camera (FOC) was built by the European 
Space Agency, but the COSTAR corrector subsequently installed aboard HST 
to correct the figuring error of its main mirror destroyed the size matching of 
the Lyot stop in the FOC coronagraph and made it inoperable. Also, the aper- 
ture size had been shrunk from 3m to 2.4m, and the telescope rotation proved 
difficult to achieve. The FOC could however have succeeded since exo-planets 
brighter than expected were recently found by radial velocity measurements 
(Mayor and Queloz 1995, Marcy and Butler 1996, Delfosse et al. 1998). 

The rotation modulation principle was later adopted, together with a form of 
coronagraphic field-selective attenuation, by Braccwell and McPhie (1979) for 
their proposed interferometer exploiting the improved contrast of a planet in 
the infrared. The concept has been refined by Leger et al. (1996) for the DAR- 
WIN proposal, and also adopted by Angel and Woolf (1997) for their Terrestrial 
Planet Finder. Both instruments are intended to observe in the 10/im infra-red. 
At this wavelength, the Earth's luminosity is quoted to be 7.10 6 times, or 17 
magnitudes, fainter than the Sun's. 

Coronagraphic imaging is now also considered for the Next Generation Space 
Telescope (Gezari et al. 1997, Moutou et al. 1998, Trauger et al. 1998, Le Fevre 
et al. 1998). The added refinement of "dark-speckle imaging" can further im- 
prove the detection sensitivity 10 to 1000 times (Boccaletti et al. 1998a), while 
the use of a phase mask and the active attenuation of residual starlight in the 
annular field push the sensitivity towards the exo-planet detection threshold at 
visible wavelengths. 

In this article we analyze and confirm the suggestion (Labeyrie 1999b) that 
these techniques of apodized imaging are also usable with large multi-aperture 
interferometric arrays operating in the densified-pupil imaging mode. Indeed, it 
was recently shown that such arrays operated in this mode can produce snap- 
shot images in a narrow field (Labeyrie 1996). We have studied the corona- 
graphic schemes applicable to such imaging conditions and find them suitable 
for exo-planet observing at infra-red wavelengths. At visible wavelengths they 
can provide higher angular resolution than NGST coronagraphy, allowing ob- 
servations of fast orbiting exo-planets. 

Our numerical simulations of imaging with the "Exo-Earth Discoverer" (EED) 
sketched in Fig. [l] verify the expected imaging performance. The EED is a 36- 
element space interferometer with size D of the order of 100m. Having free- flying 
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telescopes (Labeyrie 1985, Bcly et al. 1996), it can be considered as a possible 
precursor of much larger instruments such as the 150 km Exo-Earth Imager pro- 
posed for the longer-term goal of making exo-planet portraits (Labeyrie 1999b). 

2 Science 

The proposed instrument is intended for imaging, in the infra-red, and possibly 
the visible, faint circumstellar environments with their diffuse component and 
their point sources, brown dwarf companions or exo-planets. Brown dwarfs are 
expected to be 10 7 to 10 8 times fainter in the optical range and 10 5 times in the 
near infra-red. Extra-solar planets are 10 6 times fainter than their parent star 
at 10 microns. Their detection may also be affected by the presence of extended 
emission, both zodiacal and exo-zodiacal in their vicinity. 

Jets and disks, as well as extreme cases of circumstellar activity such as in 
SS 433 (Spencer 1979) are obviously also relevant to high-resolution coronag- 
raphy. Extended sources in the vicinity of stars are of interest, especially at 
the birth and death stages of the star life: observations of nebulosities leading 
to the formation of a planetary system, ejected envelopes, accretion disks in 
binary systems, disks of small bodies and planetary debris will all benefit from 
high-resolution imaging coronagraphy, as well as high-velocity jets and super- 
nova remnants. 

High-angular resolution can also access the inner parts of active galactic 
nuclei, inside the dust torus, and coronagraphy can improve the rejection of the 
bright unresolved core. In objects such as M81, QSO's and possibly gamma-ray 
bursters, where a central black hole is suspected, both the transverse and radial 
components of the gas velocity may become observable with milli-arcsecond 
angular resolution and the associated field spectrography. Also interesting is 
the observation of rings around supernovae in nearby galaxies. 

3 Concepts 

3.1 Pupil densification 

Fizeau interferometers, the equivalent of giant telescopes having a sparse mosaic 
mirror, produce an image, but its quality degrades catastrophically when the 
aperture size becomes much larger than the sub-apertures. A usable image can 
then be retrieved by dcnsifying the exit pupil, i.e. distorting it to increase the 
relative size of the sub-pupils. Such instruments, which may be called "hyper- 
telescopes", evade a requirement long believed to be a "golden rule of imaging 
interferometry" (Traub and Davis 1982, Beckers 1997), namely that the exit 
pupil be identical to the entrance pupil. Instead, one preserves only the ar- 
rangement of sub-pupil centers, while magnifying each sub-pupil with respect 
to the inter-pupil spacings, thus making the exit pupil more densely packed than 
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the highly diluted entrance pupil (Fig. ||). 

The result is a combined image where the interference pattern is magnified, 
with respect to the broad diffraction peak contributed by the sub-apertures. 
This broad peak acts as a window which limits the field but concentrates the 
energy in the useful central part of the interference pattern. The pupil den- 
sification mixes two different scales in the image: the low-resolution scale of 
the sub-images, or images contributed by a single sub-aperture, and the high- 
resolution scale of the interference pattern, which retains the classical convo- 
lution behaviour on extended objects. When the array is phased on a point 
source, a peak appears in the interference function, thus providing a direct im- 
age within the window. 

Depending on the baseline redundancy, snap-shot exposures obtained in this 
way can contain between N and N 2 resolved elements, if N is the number of 
apertures. Simulations of an exo-planet's image, showing clouds and continents, 
have been obtained assuming a 48 or 150- aperture "Exo-Earth Imager" , a larger 
version of the EED with 150 km size (Labeyrie 1999b). 

Here we concentrate on the shorter-term goal of obtaining unresolved im- 
ages of exo-planets, and the other goals mentioned in the previous section. It 
requires a smaller instrument, spanning 50m to 500m, so that the central source 
be unresolved, as required for phase-mask coronagraphy (Roddier and Roddier 
1997, Guyon et al. 1999) and for dark-speckle imaging (Labeyrie 1995, Boc- 
caletti et al. 1998a, Boccaletti et al. 1998b). 

For a sparse array geometry allowing full pupil densification, we assumed an 
entrance aperture shaped like an "exploded" hexagonal paving (Fig. ||). The 
case of 36 elements, with concentric rings of 6, 12 and 18, was mostly used in 
the simulations, but a wider imaging field is obtainable with more elements, at 
no loss if the collecting area is conserved. Once densified, the 36-element pupil 
resembles the aperture of the Keck telescope (see Fig. ||). 

3.2 Phase-mask coronagraphy with multi-aperture arrays 

Classical coronagraphs since Lyot (1939) have an occulting mask in the first fo- 
cal plane of a telescope, and, in a relayed pupil plane, a "Lyot stop" diaphragm. 
Slightly smaller than the pupil, it removes the double edge ring caused by diffrac- 
tion, which contains most of the light propagated from the Airy rings (Malbet 
1996). The focal mask is opaque, and has to cover at least the 4 or 5 central rings 
of the Airy pattern for efficient nulling. What remains in the final re-imaged 
field is a few attenuated Airy rings, becoming broken into random speckles at 
increasing distances from the vanished Airy peak. The speckles result mostly 
from the residual bumpiness of the telescope's mirror, but obscurations of spi- 
ders, segmentation, etc., on the pupil can contribute. 

With Roddier and Roddier's proposal (1997) of a transparent phase-shifting 
mask, smaller than the Airy peak (ps 0.5-^), similar extinction is achievable, 
but the usable field extends closer inwards, down to the very edge of the Airy 
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peak. The phase mask is basically narrow-band. It can however be multiplexed 
for a wide spectral coverage by using a polychromatic Bragg hologram as a re- 
flective phase mask (Labeyrie 1999b). 

When a star's Airy peak is properly focused and centered on a phase mask, 
most light is diffracted outside the geometric pupil, and is thus removable by a 
diaphragm, the Lyot stop, before re-imaging the focal plane. If made to match 
exactly the geometric pupil, including any central obscuration, spider arms and 
segmentation gaps, the Lyot stop can be nearly optimal. In the presence of 
obscurations or mirror segmentation such as shown on top of Fig. ||, this is 
another advantage with respect to conventional Lyot coronagraphy, where the 
bright double fringe along edges requires masking a significant part of the geo- 
metric pupil area (typically 20% of the diameter is occulted). 
Phase-mask coronagraphy for single telescopes such as the NGST (Moutou et al. 
1998) can be extended to the case of a densified-pupil interferometer. Whether 
the wavefront focused at the entrance of the coronagraph has natural continu- 
ity or is a dense mosaic of wavefront patches coming from widely spaced sub- 
apertures in the entrance aperture makes no difference indeed. A non-resolved 
star is therefore "nulled" identically in cither case. 

Another new coronagraphic scheme, the Achromatic Interfero-Coronagraph 
(Gay J. and Rabbia 1996, Rabbia et al. 1998), also offers the capability of imag- 
ing a companion very close to the field center, and uses an ingenious solution 
for achromatism. Although not included in the present simulations, it would be 
of interest to assess its performance in the EED context. 

4 Numerical simulations 
4.1 The algorithm 

We have performed numerical simulations of the proposed interferometer and 
coronagraph, using various occulting and apodizing schemes. In accordance 
with the theory of densified-pupil interferometry (Labeyrie 1996), the complex 
amplitude distribution in the exit pupil is considered as a convolution of a sub- 
pupil with a "fakir board", i.e. an array of Dirac peaks representing the centers 
of the sub-pupils (the peaks can be complex numbers to represent phasing er- 
rors) . The combined focal image of a point source is therefore a product of an 
interference function, the squared modulus of the fakir board's Fourier trans- 
form, and a diffraction function similarly transformed from the sub-aperture. 

If the array is correctly phased, the interference function resembles the clas- 
sical Airy pattern with its central peak, although the outer rings here become 
broken into speckles. The usual phasing algorithms serving for adaptive optics 
in monolithic telescopes, using wavefront slope or curvature measurements, are 
not applicable to long baselines, but other approaches using sharpness criteria 
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have been found usable on Earth and in space (Pedretti and Labeyrie 1999, 
Pedretti 1999). 

An extended incoherent source convolves the interference function, but affects 
negligibly the diffraction function if the interferometer is highly diluted. The 
interference function may therefore be considered as a spread function, while the 
diffraction function behaves like an envelope or field-limiting window applied to 
the convolved pattern. 

For simulating the phase-mask coronagraph attachment, we introduced a 
180° phase shift in the central part of the star's Airy-like interference peak. 
The mask size is quite critical and was adjusted by trial and error for minimal 
residual star light in the annular field. The Fourier transform of the result- 
ing complex amplitude distribution is then multiplied by the Lyot stop, in most 
cases made to match exactly the geometric pupil, and again Fourier transformed 
to obtain the final image of the star. 

Adding the planet's image requires a similar calculation to be made for the 
planet. In the plane of the phase mask, this involves the same interference 
function, but attenuated and shifted off-axis before multiplication with the un- 
shifted diffraction function and the phase mask function. Since the planet's Airy 
peak falls outside of the phase mask, the mask affects very little the planet's 
final image. 

Various sizes and shapes of the mask and pupil diaphragm were also experi- 
mented for optimizing the nulling of the star and its feet. 

4.2 Comparison of various pupil configurations 

We have considered diluted arrays of 36 elements, having either circular sub- 
apertures or, better, hexagonal ones allowing a full densification in the exit 
pupil. Although the interference functions are identical in both cases, we find 
that a darker coronagraphic field is achievable with the fully densified pupil. 
We therefore adopted this optimal pupil configuration for simulating the coro- 
nagraph in the following. The pupils and corresponding images are shown in 
Fig.§ 

When working with hexagonal sub-apertures and full densification, the first- 
order secondary peaks of the star's interference function are exactly located on 
the envelope's dark ring. The field area, expressed on the sky in units of the 
array's Airy area (squared Airy radius) is approximately equal to the number 
of apertures. The secondary peaks of the star's interference function have the 
same intensity as the central peak, but have first-order lateral chromatism. Any 
stellar companion or planet present in the main field or in the extended field 
also has side-peaks. In this respect, in monochromatric light, the central field 
always contains a replica of the planet's image if the central image is outside 
(see Fig. 5a). 
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4.3 Phase-mask size and shape 

Here the sub-apertures of the interferometer are assumed perfectly phased. 
The maximal nulling of starlight by the phase mask is reached when the ampli- 
tudes with opposite phases are balanced in the pupil, and this requires a careful 
adjustment of the phase mask diameter (Fig. ||, top). A circular phase mask 
appears better than an hexagonal one. The calculated intensity profile in the 
coronagraphic image of the unresolved star (Fig. ||, bottom) shows that the av- 
erage attenuation of the annular held achieved by the phase-mask coronagraph 
is of the order of 20 to 50, a value comparable to that obtained by Roddier 
and Roddier (1997) for a monolithic telescope. Although not a large gain, this 
brings the average halo level to w 2.10~ 4 relative to the peak of the unmasked 
star image, a welcome improvement before applying the further darkening steps 
discussed below. 

Residual phase errors, both the intra-aperture errors caused by the bumpiness 
of the mirror elements and the inter-aperture errors caused by imperfect ad- 
justment of the optical path differences among them, have the same effect as in 
monolithic telescopes: the feet of the diffraction pattern become distorted and 
intensified, both before and after the coronagraph. The intra-aperture bumpi- 
ness diffracts mostly outside of the imaging held and therefore contributes little 
to the image degradation. In the 10 micron infra-red, the inter-aperture path 
differences can be kept small compared to the wavelength. Our simulations 
indicate that wavefront bumpiness contributes little to the background. If its 
scale size is 10 times smaller than a sub-aperture and its RMS amplitude 59nm, 
corresponding to A/8 at 0.5/zm and A/170 at 10/xm), the bumpiness contribu- 
tion to the 2.10~ 4 background is only 10~ 8 . Subtracting a reference star image, 
obtained with a different bumpiness, therefore leaves a speckled background 
at 10~ 8 level (Fig. [|, bottom), and it can be smoothed if the actuators are 
re-ajusted many times during a long exposure. In the dark-speckle mode, the 
analysis is performed differently. 

4.4 Effect of full or partial pupil densification 

Another critical parameter is the amount of densification in the exit pupil, and 
numerical simulations have again served to characterize its effect (Fig. [|). We 
varied the sub-pupil size and calculated the nulling in the annular coronagraphic 
field, extending from the second dark ring of the masked Airy peak p mi outward 
to the first dark ring of the sub-aperture's Airy peak pf. With 36 elements, its 
area is 25 times larger than the area of the first dark ring p u in the unmasked 
Airy peak (Fig. ||, bottom). The residual intensity reaches 2.10~ 4 with full 
densification. It increases 10 times for 50% density. The optimal phase-mask 
size depends on the density and may be adjusted for optimal nulling 
Narrow gaps between elements, amounting to a few percents of the element 
size, are difficult to avoid in practice, but will not degrade significantly the 
field darkness. This is another advantage with respect to the classical Lyot 
coronagraph, where such amplitude patterns on the pupil can be disastrous. 
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5 Application to exoplanet detection 



5.1 Choice of wavelength 

The planet/star contrast improves a lot from the visible to the 10/xm infra-red, 
owing to the drop of the star's Planck function and the peaking thermal emis- 
sivity of the planet. Also, residual wavefront bumpiness causes smaller phase 
shifts and less scattered light at longer wavelengths. However, the angular sepa- 
ration of an 8m space telescope such as the NGST becomes sufficient to separate 
planets from their star in the visible, although not in all cases, and provides in 
principle the same planet detection performance as an EED of identical collect- 
ing area. 



5.2 Performance of direct imaging 

Optimal data reduction would involve a combination of the dark hole and dark 
speckle techniques mentioned below. The simulation routines are demanding 
of computer time however, and we have here simulated long-exposure images, 
a less sensitive method which smoothes the boiling speckles and thus cannot 
exploit the nulling achieved in the dark speckles, but nevertheless demonstrates 
the detectability of planets. The results shown in Fig. || assume 36 phased 
apertures of 0.6m, for a total collecting area similar to that of DARWIN (about 
40 m 2 ) (Leger et al. 1996, Mennesson and Mariotti 1997). We have simulated 
a twin of the solar system located at 20pc, with Venus, Earth and Mars near a 
m v = 6.33 star. At 10/^m the stellar magnitude is 11.5. The detector noise and 
photon noise were also included in the simulated 10 hours exposure. Planets 
w 7.10 6 times fainter than the parent star are directly detected after subtracting 
a reference image. 

Phasing errors, modelled to represent the bumpiness of the mirror segments, cor- 
responding to A/170 rms, are included both for the planetary system and the 
reference star. A long exposure was fabricated by adding 60 images obtained 
with independant random phase maps, and this was repeated with different 
phase maps for the reference star. With larger wavefront errors of A/80 (at 
10/im), the scattered light level is about 10 -7 and is therefore higher than the 
Airy peak of Mars which remains undetected. 

We have carried out the simulation at 3 different wavelengths (9.5, 10 and 10.5 
/im). Secondary peaks feature a I s * order chromatic dispersion and cannot be 
confused with the white primary peaks. Due to the interferometer geometry, an 
off-axis source has 3 peaks at 120° (1 primary and 2 secondary). 
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5.3 Further gain with dark-hole and dark-speckle tech- 
niques 

Malbet et al. 1995, and Trauger et al. 1999, have shown how minor corrections 
of the wavefront shape, applied with active optics, can further improve the dark- 
ening of a coronagraphic image. Using a coronagraphic exposure, the dark-hole 
algorithm derives a phase correction map. Once applied to the actuators, it 
improves the average darkness of the selected annular field. 
The pattern of residual speckles changes randomly after each iteration, a conse- 
quence of noise in the servo loop. Although these speckles are here much fainter 
and slower than those caused by residual turbulence on Earth, the situation is 
similar, and dark-speckle imaging (Labeyrie 1995, Boccaletti et al. 1998a) can 
similarly improve the star-light rejection if thousands of short exposures can be 
exploited. The actuators can be re-adjusted every 10s for example, and frozen 
in-between during the exposures. 

Dark-speckle observations performed on ground-based single-aperture telescopes 
(Boccaletti et al. 1998b) with a CP20+ photon counting detector (Abe et al. 
1998) have verified the theoretical expectations. Their extrapolation to the 
situation considered here suggests that exoplanets (Am(10/xm) w 16 to 19, 
Am(lfim) w 22 to 25) can be detected with exposures lasting hours to a few 
tens of hours. A limitation is the imperfect cleaning of the image, with a few 
stellar speckles survfiving near the center. These can be removed with a refer- 
ence star, but perhaps also by modifying the dark-hole algorithm. 

5.4 Signal to noise ratio 

In addition to the residual pattern of star-light, the image provided by a densified- 
pupil interferometer, contains light from the zodiacal and exo-zodiacal clouds. 
The part of these extended sources which lies outside of the narrow imaging 
field, but within the resolution patch of the sub-apertures, contaminates the 
image through the feet of the interference function. 

An exo-planet peak in the image, with J p photo-events detected every sec- 
ond, is thus contaminated with : 

1- J s events/s from the residual speckled halo of starlight; 

2- J ez events/s from the image of the exo-zodiacal cloud; 

3- J' ez events/s from the feet of the exo-zodiacal image (the cloud being larger 
than the field, there is a contamination of the field from the missing parts of 
the cloud, through the feet of the spread function); 

4- J z events/s from the zodiacal cloud; 

5- Jt events/s from the telescope thermal emission. 

To calculate these quantities, we consider an entrance aperture of size D 
with N sub-apertures of size d. It receives I p photons/s per square meter from 
the planet; I s from the star; L z from the zodiacal cloud, per steradian; and L ez 
per steradian from the exo-zodiacal cloud of angular size <j), assumed unresolved 
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by a sub-aperture. 

Since the planet's peak concentrates most entering photons, J v — I p Nd . 
With the long-exposure mode using a reference star, the subtracted stellar con- 
tribution at the planet peak may be approximated by J s = I s Nd 2 g. 
The radial profile g = p~ a of starlight is estimated if p is the axial distance and 
a the log slope of the coronagraphic halo, which is of the order of a w 1.5 ~ 2 
according to the simulation results of Fig. |^ (bottom), a decreases slowly as 
the number of speckle patterns is increased, as does the speckle noise. 
The remaining components may be expressed as : 

J z = L Z X 2 in the zodiacal cloud ; 

J ez = L ez \ 2 Nd 2 D~~ 2 in the image of the exo-zodiacal cloud; 
J' ez = L ez <p 2 d 2 in the halo of the exo-zodiacal image; 

Jt = LtX 2 since the telescope thermal emission produces a background similar 
to the zodiacal cloud. 

The last four quantities describe the respective contributions of the central 
peak and the infinitely wide pattern of secondary peaks in the interference 
function. This pattern adds a halo to the image of the exo-zodiacal disc. A 
similar halo is also added to the image of the zodiacal continuum, and dominates 
it completely. 

These quantities are invariant when N varies at constant collecting area Nd 2 , 
except J' ez which decreases for increasing N. Increasing N, which is of interest to 
enlarge the field window, thus also improves the planet's discrimination against 
the halo of the exo-zodiacal image. In the N — > +00 limit, the sparse array 
with densified exit pupil provides the same image as a filled telescope of similar 
diameter, although attenuated in the ratio of collecting areas. 
In the long-exposure observing mode, noise components are: 

1 . photon noise P n from the added photon contributions in the same speckle 
as the planet peak, amounting to P n = (2 J s + J p + J z + J ez + J'^T , 
where T is the total exposure time (the factor 2 in the stellar term arises 
from the subtraction of a reference image) ; 

2. speckle noise S n from the stellar contribution, S n = y / 2J 2 T 2 /n s (Angel 
1994), if n s is the number of speckle patterns added in the long exposure. 
For n s > 20 the speckle noise becomes negligible; 

3. read-out noise, R n = y / 2a ron n, where a ron is the detector's read-out noise 
per speckle area, whith n exposures made on both the object and reference 
star; 

4. thermal emission within the telescope: Diner et al. (1991) calculated that 
a telescope at 70K provides, at 10/im, a 10% thermal contribution in the 
image of the zodiacal cloud. At 40K, the mirror temperature considered 
for DARWIN and TPF, the ratio is in the range 10~ 2 - 10~ 3 . Here, 
this ratio is unaffected when the Fizeau images goes through the densificr 
optics, provided its temperature is not higher than the primary mirror's. 
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The image contamination from thermal emission is therefore negligible for 
planet searching at 10 fj,m. 



The long-exposure signal/noise ratio is therefore : 



S/Need 



I p Nd 2 T 



(1) 



or 



S/Need 



I p Nd 2 T[(2J s + J P + J Z + J ez 
+ J' ez )T + 2J 2 T 2 /n s + 2a ron n}- 1 / 2 



(2) 



The primary peak from the planet becomes attenuated if it moves towards 
the edge of the field. However this causes two secondary peaks to move inwards 
and become brighter. Because their location is known a priori, they can be 
added to the primary peak. The total signal thus obtained is invariant with 
respect to planet position, the energy being conserved. And the added contri- 
butions from uniform sky background within the planet's peaks are similarly 
invariant for the same reason. The same applies to the added contributions 
from the star's residual halo if it has the same radial attenuation as the en- 
velope. The chromatic elongation of the secondary peaks does not necessarily 
affect the addition of the planet peaks since it is removable with a field spec- 
troscopy arrangement. Therefore, the S/N is little affected by a planet's motion 
in the field. 

In the dark-speckle mode, the stellar contribution to photon noise and speckle 
noise is decreased, as calculated by Boccaletti et al. (1998a). Also, no reference 
star is needed if the residual fixed speckles can be discriminated from planets by 
refining the dark-hole routine. Whether sensitive wave analyzer schemes, such 
as Zernike's phase contrast, can help remains to be investigated. 

5.5 Comparison with Bracewell nulling 

In a multi-aperture Bracewell interferometer, light from a field extending over 
the resolution patch of the sub-apertures is transmitted to the single detector 
pixel, although with a field-dependant attenuation defined by the high-resolution 
transmission map of the beam-splitter arrangement. This map has a central 
dark minimum, designed to attenuate the star, surrounded by brighter speck- 
les which can be rotated or "boiled" for generating a planet signal. Like the 
procedure of Bracewell and Mac Phie (1979), this is a single-pixel version of 
the rotation/subtraction algorithm originally proposed for exo-planet searching 
with the Hubble Space Telescope by Bonneau et al. (1975) (see also Davies 
1980). The angular transmission map is such that the average transmission r 
of the planet's light and zodiacal light is of the order of r = 0.2. 
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With respect to the J terms previously listed, the zodiacal and exo-zodiacal 
contributions must now be multiplied by Nt since the detector pixel now re- 
ceives all the energy from the part of these sources appearing in the sub- 
aperture's resolved patch. Also, the exo-zodiacal image vanishes, as well as 
the speckle noise. Thus, the signal/noise expression becomes: 

S/N BraceweH =I p Nd 2 TT{(J s + J p t + J z Nt 
+ j' ez NT)T + a^n]- 1 ' 2 

where g now describes the starlight rejection factor, expected to reach 10~ 6 in 
DARWIN. 

Increasing the number of sub-apertures N at constant collecting area NcP 
changes only the zodiacal and exo-zodiacal signals, which increase in proportion 
to N. Unlike the EED, Bracewell interferometers should thus remain limited to 
few apertures. 

Using these expressions, we compared the performance of an EED having 
36 apertures of 0.6m, and a DARWIN of equivalent area, both located at 1 
A.U. from the Sun. The instrument's transmission combined with the quantum 
efficiency of the detector is 50% at 10/xm. An Earth at 20 parsecs provides 4 
ph/m 2 in one hour of exposure. 

At lO^im, the luminosity of the solar zodiacal cloud varies between 
7,6.10- 7 W/(m 2 .sr .fj,m), at the ecliptic pole, the value used by Mennes- 
son et al. (1997) and 1, 1.10^ 4: W/(m 2 .sr.fim). We used the median value 
(7, 6. 10~ 6 W / (m 2 .sr. /im) corresponding to 21 magnitudes per square arc- 
sec. The equivalent angular size (f> of the exo-zodiacal cloud is taken as 
4> = 0.25arcsec, corresponding to 5 A.U. in the solar system, with L ez = 10L Z . 

A 30 hr observation in the narrow spectral band (A = 10/ito, AA = 0.5/xm) 
provides a 0.5 o detection with DARWIN, and a 5.8 a detection with the 
EED in the long-exposure mode with "boiling" speckles. With an achromatic 
phase-mask allowing a wider 6 — 18/xm band, the DARWIN spectroscopy band, 
the signal/noise ratio is 50.8 for the EED and 4.7 for DARWIN. 

The gain in the EED case results from two effects: 

1- the fully constructive interference forming the planet peak in the image, as 
opposed to the 20% average intensity in the transmission map of DARWIN. 

2- the fact that the planet peak can be isolated from most of the zodiacal 
contamination in the high-resolution image. 

For a given collecting area, the array size does not affect the signal/noise 
ratio. With the EED, the image is magnified, within its fixed window, when 
increasing the array size, and this can bring the planets optimally located in 
the field. Increasing the sub-aperture count N from 36 to 100, their size being 
kept constant, increases EED's signal/noise ratio nearly proportionally from 
5.8 to 16.0. Instead, DARWIN's value of 0.5 in 30 hours increases to 0.8, as the 
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square root of the aperture number. 



6 Conclusion and future work 

The prospect for coronagraphic imaging with large interferometric arrays, at 
infrared and visible wavelengths, opens vast areas of science. The simulations 
presented show that searches for exoplanets at 10/xm should be feasible once the 
formation flight of mirror elements in space becomes mastered. We recall the 
main conclusions : 

• Pupil densification allows direct imaging and coronagraphic nulling. 

• A phase-mask, replacing the classical Lyot opaque mask, cleans the image 
closer to the attenuated star, and is more tolerant to the unavoidable nar- 
row gaps between sub-pupils in the densified pupil. Both achromatisation 
of the mask and low-resolution spectroscopy (A/AA = 20) are desirable 
for analysing the images obtained. 

• At given collecting area, more sub-apertures is better in terms of image 
field and signal/noise ratio. 

• The optimal wavelength is a trade-off between resolution, stellar spill- 
off, planet contrast, background emission and detector performance. At 
visible wavelengths, a monolithic 8m telescope often resolves the planet 
from its star and then has the same planet-searching efficiency as an EED 
of identical area. 

• In addition to planet searches, direct imaging with the EED is also of 
interest for the other science goals briefly described above. 

• With respect to Bracewell's detection scheme, and its DARWIN and TPF 
proposed implementations, the image formation in the EED improves 
markedly, at equal collecting area, the discrimination of an exo-planet 
amidst the zodiacal and exo-zodiacal nebulosity. The EED however re- 
quires 36 apertures at least, but these are smaller and can be implemented 
with identical mirror segments compatible with the nano-satellite philos- 
ophy . 

The coronagraphic image is compatible with field spectrography. 

Additional simulations are needed to compare the possible exposure 
routines, assess the desirability of a fine wavefront sensor, and refine the noise 
models. More work is also needed to design optimal optical trains, including 
the coronagraphic device and its spectro-imaging attachment. 

Acknowledgments : Jean-Marie Mariotti his constructive was among the 
first to understand our ideas and on densified-pupil imaging encouraged us to 
explore their application. His premature death deprived us of his insight and 
critical appreciation, particularly regarding these results. 
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Figure Captions 



Figure U Possible optical train for the Exo-Earth Discoverer. Mirror 
elements Ml focus starlight in the focal combiner C (detail below), on field 
mirrrors FM, which form a densified image of the pupil on mirror segments 
PM. A combined image is focused in F at the entrance of a coronagraph Co. 
Mirrors PM are carried by actuators for tip-tilt and piston corrections. 

Figure |^: Interferometer apertures (top row) with disk (left) or hexagon 
(right) subpupils, and densified version of them obtained in the exit pupil 
(second row). The 100 to 200m Exo-Earth-Discoverer has a sparse entrance 
aperture, but a fully densified exit pupil such as shown at right in the second 
row. The star image (third row), with its diffractive envelope, has most energy 
in the central Airy peak, as seen in the intensity profiles (bottom). 

Figure |^: Top: Central darkening versus phase-mask size with a fully 
densified pupil of 36 adjacent hexagons, measured at the maximum ( dotted 
line) and across the area (solid line) of the Airy peak . The optimal mask size 
found is (0.555 ± 0.005)A/D. Bottom: Azimuthally averaged and normalized 
profiles of the unmasked Airy pattern (solid line), coronagraphic pattern 
(dashed line) and residual scattered pattern after subtraction of the reference 
frame (dotted line), including photon and readout noise. Both long exposures 
incorporated 20 different maps of (A/ 170 rms) wavefront bumpiness. The 
radii of the first dark ring in the unmasked (p u ), the masked (p m ) and the 
sub-aperture (pf) images are indicated. Also indicated is the level of Venus 
(dash-dot), Earth (dash-dot-dot-dot) and Mars (long dashes) multiplied by the 
sub-aperture's diffraction function . 

Figure ^: Coronagraphic attenuation of star light versus exit-pupil density. 
The pupil density is measured as the relative area of the sub- pupils in the 
exit pupil, while the attenuation is measured in the annular field (12X/D) and 
normalized to the intensity of the unmasked peak. The dotted curve indicates 
that the mask size was re-optimized for each density value, while the solid 
curve indicates that the initial value (0.555A/D) was maintained . Numbers 
indicate the sizes of the re-optimized masks in units of X/D. 

Figure ^| Coronagraphic image of a solar system twin observed at 20pc 
(m v — 6,33) with the Sun attenuated to detect Venus, Earth and Mars. The 
magnitude differences at 10 pm. are respectively 16.4, 17.1, 19.8 relative to 
the parent star. The simulated array has 36 telescopes of 0.6m, providing a 
collecting area similar to a DARWIN (40m 2 ). To discriminate the dispersed 
secondary peaks from the white primary peaks, the final image was obtained 
with a combination of 3 wavelengths (9.5, 10 and 10.5/im). 
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In (a), the star is removed to show the Venus, Earth and Mars components of 
the image, with their primary (V, E, M) and secondary (V\, V2, E\, E 2 , Mi, 
M2) peaks. The peripheral attenuation by the diffraction envelope makes the 
secondary peaks of Mars brighter than the primary peak, since they are closer 
to the field center. No photon noise is present here. 

The separations are 1.75A/D, 2A3X/D and 3.69A/Z?, corresponding, for a 
100m array, to 36, 50 and 76 milli-arcsecond at 10 /im, respectively. A circle 
indicates the field limit, at the first dark ring of the diffraction envelope of a 
single 0.6m-aperture. On the sky, the corresponding field diameter is 12X/D, 
amounting to 0.24 arc-second at 10 (im with a 100m array of 36 elements. 
Image (b) shows a 10 hours exposure with subtraction of a reference star 
image similarly obtained. Photon noise and readout noise (100 e-/pix/frame) 
were included. A perfect wavefront was assumed. Image (c) is similar, but 
with 59 nm RMS wavefront bumpiness (A/170), of lOmn lifetime. Two 
different sets of 60 phase maps served for the object and the reference 
star. This level of bumpiness has little effect in the usable field, within the 
white circle. A further gain in sensitivity is achievable in the dark-speckle mode. 



16 



References 



[I] Abe L., F. Vakili, I. Percheron, S. Hamma, J. P. Ragey and A. Blazit 1998. 
in Proc Catching the perfect wave. Albuquerque June 1998 

[2] Angel J.R.P. and N.J Woolf 1997. An imaging nulling interferometer to study 
extrasolar planets. Astrophys. J. 475, 373-379 

[3] Beckers J. 1997. In Instrumentation for Large Telescopes. (JJ Rodriguez 
Espinosa Ed.) Cambridge 

[4] Bely P.-Y., R.J. Laurance, S. Volonte et al. 1996. Kilometric baseline space 
interferometry. In SPIE Proc. Space Telescopes and Instruments IV (P.-Y. 
Bely, J.B. Breckinridge Eds.), pp. 59-73. Vol 2807 

[5] Boccaletti A., A. Labeyrie and R. Ragazzoni 1998. Preliminary results of 
dark-speckle stellar coronagraphy. Astron. Astrophys. 338, 106-110 

[6] Boccaletti A., C. Moutou , A. Labeyrie , D. Kohlcr and F. Vakili 1998. 
Present performance of the dark-speckle coronagraph Astron. Astrophys. 
Suppl. Ser. 133, 395-402 

[7] Bonneau D., M. Josse and A. Labeyrie 1975. In proc. Utrecht Symp. Image 
processing techniques in Astronomy (de Jager/Nieuwenhuijzen eds) pp. 403- 
409. D. Reidel, Holland 

[8] Bracewell R. and R.H. Mac Phic 1979. Searching for nonsolar planets. Icarus 
38, 136-147 

[9] Davies D.W. 1980. Direct Imaging of Planetary Systems around Nearby 
Stars. Icarus 42, 145-148 

[10] Delfosse X., T. Forveille, M. Mayor, C.Perrier, D. Naef and D. Qucloz 1998. 
The Closest extrasolar planet. A giant planet around the M4 dwarf G1876. 
Astron. Astrophys. 338, L67-L70 

[II] Diner D.J., E.F. Tubbs, S.L. Gaiser and R.P. Korcchoff 1991. Infrared 
imaging of extrasolar planets. J. Brit. Interplanetary Soc. 44, 505-512 

[12] Gay J. and Y. Rabbia 1996. Principe d'un coronographe interferentiel. CR. 
Acad. Sci. Paris t. 332, Scric II b, 265-271 

[13] Gezari D. Y., A. Crotts, W. Danchi ct al. 1997. An optimized corona- 
graph for Next Generation Space Telescope. Columbia University and Ball 
Aerospace Systems Division, proposal in response to NRA 98-GSFC-l 

[14] Guyon O., C. Roddier, J.E. Graves, F. Roddicr, S. Cuevas and C. Espejo 
1999. Nulling coronagraph. Submitted to Applied Optics. 

[15] Labeyrie A., B. Authier, T. de Graauw, E. Kibblewhite and G. Weigelt 
1985. TRIO: A kilometric array stabilized by solar sails. In ESA Colloq. Kilo- 
metric Optical Arrays in Space pp. 27-33. Vol 226 



17 



[16] Labeyrie A. 1995. Images of exo-planets obtainable from dark speckles in 
adaptive telescopes. Astron. Astrophys. 298, 544-548 

[17] Labeyrie A. 1996. Resolved imaging of extra-solar planets with future 10- 
100km optical interferometric arrays. Astron. Astrophys. Suppl. Ser. 118, 
517-524 

[18] Labeyrie A. 1999, Direct Searches: Imaging, Dark-Speckle and Coronagra- 
phy. in NATO ASI Cargese summer-school Planets outside the Solar System 
pp. 261-279. (Mariotti J.-M. and Alloin D. Eds.) 

[19] Labeyrie A. 1999. Exo-Earth Imager for exoplanet snapshots with resolved 
details, to appear in ASP Conf Series Working on the Fringe. (S. Unwin and 
R. Stachnik, Eds.) Dana Point, CA, May 24-27, 1999 

[20] Le Fevre O. 1998. Study of Payload Suite and Telescope for the NGST. in 
response to ESA AO/ 1-3405 

[21] Leger A., J-M. Mariotti, B. Menesson, M. Ollivier, J.-L. Pugct, D. Rouan 
and J. Schneider 1996. Could we search for primitive life on extrasolar planets 
in the near future? the DARWIN project. Icarus 123, 249-255 

[22] Lyot B. 1939. A study of the solar corona and prominences without eclipses. 
Mon. Not. Roy. Astron. Soc. 99, 580-594 

[23] Malbet F., J.W. Yu, and M. Shao 1995. High-Dynamic-Range Imaging 
Using a Deformable Mirror for Space Coronography. Publ. Astron. Soc. Pac. 
107, 386-398 

[24] Malbet F. 1996. High angular resolution coronography for adaptive optics. 
Astron. Astrophys. Suppl. Ser. 115, 161-174 

[25] Marcy G.W. and R.P. Butler 1996. A Planetary Companion to 70 Virginis. 
Astrophys. J. 464 L147-151 

[26] Mayor M. and D. Queloz 1995. A Jupiter masss companion to a solar-type 
star. Nature 378, 355-359 

[27] Mennesson B. and J.-M. Mariotti 1997. Array configurations for a space in- 
frared nulling interferometer dedicated to the search for Earth-like extrasolar 
planets. Icarus 128, 202-212 

[28] Moutou O, A. Boccaletti and A. Labeyrie 1998. A Coronagraphic Dark- 
Speckle Imager for the NGST. In 34th Liege Astrophysics Colloq. NGST: 
Science Drivers and Technological Challenges pp. 211-216. 

[29] Pedretti E. and A. Labeyrie 1999. A hierarchical phasing algorithm for 
multi-element optical interferometers, Astron. Astrophys. Suppl. Ser., 137, 
1-9 



18 



[30] Pedretti E. 1999. Phasing algorithms for multi-element optical interferom- 
eters, to appear in ASP Conf Series Working on the Fringe. (S. Unwin and 
R. Stachnik, Eds.) Dana Point, CA, May 24-27, 1999 

[31] Rabbia Y., P. Baudoz and J. Gay 1998. Achromatic Interfero Coronagraphy 
and NGST. In 34th Liege Astrophysics Colloquium pp. 279- NGST: Science 
Drivers and Technological Challenges pp. 279-284 

[32] Roddier F. and C. Roddier 1997. Stellar coronograph with phase-mask. 
Publ. Astron. Soc. Pac. 109, 815-820 

[33] Spencer R.E. 1979. A radio jet in SS433. Nature 282, 483-484 

[34] Traub W.A. and W.F. Davis 1982. Coherent optical system of modular 
imaging collectors (COSMIC) telescope array - Astronomical goals and pre- 
liminary image reconstruction results. In SPIE Proc. International Conference 
on Advanced Technology Optical Telescopes pp. 164-175. Vol 332 

[35] Trauger J., D.E. Backman, C.A. Beichman et al. 1998. Strategies for High- 
Contrast Imaging with NGST. proposal in response to NRA 98-GSFC-l 



19 




Figure 1: 



20 




Figure 2: 



21 





5 10 

Radial distance in A\D units 



Figure 3: 



22 




Figure 4: 



Figure 5: 



23 



